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The vapor pressure of water at its triple point was measured with exceptionally high accuracy by 
realizing it with a special apparatus and measuring the pressure with the NBS precision mercury mano- 
meter. The vapor pressure apparatus had a system for circulating the liquid water. Actual triple point 
conditions were established with a thin sheet of freshly distilled liquid flowing down over an exposed 
mantle of ice frozen on a vertical well. This technique reduced non-volatile contaminants and the 
vapor was repeatedly pumped to remove accumulated volatile contaminants. A diaphragm pressure 
transducer was used to separate the water vapor from the helium used to transmit the pressure to the 
manometer. The value found for the vapor pressure of water at its triple point was 611.657 Pa with an 
uncertainty of ± 0.010 Pa from random errors, computed at 99 percent confidence limits. The systematic 
errors are estimated to be insignificant relative to the random errors. 
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1 . Introduction 

Water is the stuff of life— the working medium for 
power generation— the great highway of commerce — 
the determining factor in the weather. It has such 
importance to mankind that probably no other single 
compound has received more attention from the scien- 
tific world. Its properties have been measured and re- 
measured; thermodynamic tables have been con- 
structed and repeatedly revised, and further work of 
this type is even now commanding worldwide atten- 
tion. In part, this kind of activity continues because 
water is not an easy substance to deal with. It is almost 
a universal solvent, often reacts with the materials 
of an experiment, and is susceptible to significant 
changes in its isotopic composition. Consequently, 
high accuracy in the measurement of many of the prop- 
erties of water is hard to achieve, and more effort is 
needed to extend or improve our knowledge of them. 

It is our purpose to measure with the highest 
possible accuracy the vapor pressure of water for 
— 5 °C < *68 < 100 °C. We chose to concentrate ini- 
tially on the vapor pressure of the triple point, which 
was known with poor accuracy. Although a "triple 
point" is a precise designation for a phase diagram, it 
is not likely to be an actual point in its physical realiza- 
tion. Thermodynamically, a triple point is a system 
containing matter in three phases with each component 
in equilibrium throughout the system, i.e., the chemical 
potential of each component is the same in every phase. 
If all three phases should happen to be in contact 



simultaneously (which is not a requirement of the 
definition), the contact will be along a line producing a 
closed curve, reducing to a point only in the limit. 
Furthermore, there are several triple points of water. 
We shall use the unqualified expression, however, to 
denote the equality of the chemical potential of pure 
water in a system containing the three phases: ice I, 
liquid and vapor. This is the only triple point of water 
involving the 3 states of matter and the one of most 
importance. Not only is it the terminus of the ice I- 
vapor pressure curve and the beginning of the equilib- 
rium liquid-vapor pressure curve, 1 but also its tempera- 
ture, which can be expected to be invariant by the 
phase rule, has been assigned a value that, in principle, 
establishes all thermodynamic temperatures. This 
unique situation makes a highly accurate value of its 
vapor pressure, which is also invariant, especially 
interesting, because it can provide the foundation for 
exceptionally accurate thermodynamic calculations 
for water. This is true for this temperature region in 
particular because thermodynamic temperatures are 
inherently known with their highest accuracy in the 
region of the defining value. 

Most values reported for the vapor pressure of the 
triple point have been derived by interpolation or 
extrapolation of measurements of the two-phase 
equilibria between the liquid and vapor as a function 



1 Also it is the terminus of the ice I-liquid curve, but this is of no interest for vapor pres- 
sure measurements. 
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of temperature. The most significant of these are the 
values of Thiesen and Scheel, 610.90 Pa [l], 2 Scheel 
and Heuse, 610.90 Pa [2] and 610.55 Pa [3] and of 
Douslin, 612.1 Pa [4]. The average of the earlier results 
fitted along with data over a large temperature range 
is perhaps best represented, in the sense of a con- 
sensus, by the accepted value of the triple point in 
present day steam tables, viz., 611.2 Pa. 

For equal accuracy of pressure measurements, 
those workers who actually realized the triple point 
probably attained greater accuracy in the measured 
vapor pressure. There are only two papers reporting 
such measurements: 

1. The values of Prytz, published in 1931 [5], had 
an average of 4.5867 mm of Hg pressure, or 611.51 
Pa. 3 He realized the triple point statically, i.e., an 
incomplete layer of ice was formed on the surface of a 
quiescent pool of water. The pressure was measured 
by a mercury manometer using interference techniques 
to locate the menisci with 0.5 /xm Hg (0.065 Pa) 
imprecision, but the values of the vapor pressure 
varied by as much as 1.7 ptm Hg (0.22 Pa). Prytz was 
of the opinion that the average value was about 2 /xm 
Hg (0.26 Pa) too high. 

2. L. Besley and G. Bottomley published two values, 
611.29 Pa obtained from direct measurement and 
611.11 Pa derived by interpolation [6]. A special 
mercury manometer was used for the pressure meas- 
urements, with the height of the mercury column 
determined cathetometrically. The vapor pressure 
cell was connected directly with the manometer, as 
was also the case in the measurements by Prytz. 
Considerable effort was made to eliminate dissolved 
air in the filling of the cell. It was operated statically, 
being carefully thermostated to realize the 3-phase 
equilibrium at the surface of the liquid. The total 
uncertainty of either value was not given. The value 
derived from direct measurements was based on 148 
determinations which varied between 4.579 and 4.590 
Torr, and one might thereby regard the variation of 
±5.5 m Torr as equivalent to the imprecision at the 
99.4 percent confidence level. The residual standard 
deviation for the polynomial, from which the inter- 
polated value (4.5837 Torr) was derived, was 1.7 m 
Torr. The two imprecisions are not inconsistent, 
insofar as three residual standard deviations would 
approximate the imprecision of the interpolated value 
at the 99+ percent confidence level, except that in 
addition the uncertainty would be greater because the 
triple point was near the lower extreme of the range of 
the fit. 

Our own method also realized the triple point, but 
in a dynamic system. Liquid water is distilled onto a 
mantle of ice and falls into a well to be pumped back 
to the still. This method helps to continuously reduce 
the concentration of dissolved gases which come out 
of solution in the still and in the pressure cell and are 
pumped off periodically. Since the evaporation from 



the still is quiescent, nonvolatile impurities are retained 
in the boiler. To maintain the purity of the water vapor 
in the triple point cell, and of the helium in the mano- 
meter line, the vapor pressure system is separated 
from the manometer line by a diaphragm pressure 
transducer. 

2. Equipment 

Improved accuracy of the measurements of the vapor 
pressure of water at the triple point can only be 
attained by: 

(1) Making more accurate pressure measurements. 

(2) Realizing the triple point better. 

(3) Improving the purity of the water, especially 
with regard to gases in solution. 

The measurements were undertaken in the NBS 
Gas Thermometry Laboratory because the first problem 




2 Figures in brackets indicate the literature reference at the end of this paper. 

3 Based on the equivalence of 760 mm Hg= 101325 Pa. 



FIGURE 1. Schematic of equipment. 
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was essentially solved by using the NBS precision 
mercury manometer to measure the pressure. The 
second and third problems were dealt with by the con- 
struction of an elaboration of the triple point cell, 
and by the use of a metal diaphragm pressure trans- 
ducer. The arrangement of the equipment is shown 
schematically in figure 1, where M is the manometer, 
D is the diaphragm, and the remainder is the varying 
and the triple point cell. We shall discuss the equip- 
ment in terms of these three main elements. 

2. 1 . The Manometer 

The principles and many of the details of the NBS 
precision mercury manometer are given in an earlier 
paper [7]. Essentially it is a W-tube manometer, with 
the locations of the crowns of very large menisci pre- 
cisely reproduced by the use of capacitance measure- 
ments, and with the height of the mercury column 
accurately measured by end length standards (gage 
blocks). The instrumentation of the manometer was 
designed so that a null signal on a capacitance bridge 
was observed when the gas pressure in the lower cells 
balanced the pressure from the mercury column plus 
the vapor pressure of the mercury in the upper cell. 
The maintenance of the pressure was made automatic 
by using the output of the bridge to control a heater in 
a large ballast (a "thermal injector") so as to restore 
the pressure. 

The pressure to be measured is between 611 and 612 
Pa, corresponding to a height of the mercury column 
close to 4.6 mm. We bought a chromium carbide gage 
block of that length in the highest quality offered. It 
was carefully calibrated by J. S. Beers, Deputy Chief 
for Length in the Dimensional Technology Section. 
In figure 2 the interference pattern of this block indi- 
cates very good flatness and parallelism, with one 
corner slightly low relative to the gaging point. These 
properties make meaningful the high precision of the 
calibration measurement, which had an estimated 
standard deviation of the mean of 2.6 X 10 ~ 6 mm. This 
includes the variability due to wringing. 

Because it is impracticable to obtain a gage block of 
exactly the correct length (particularly in advance of 
knowing its value) some modification in our usual 
method of measuring pressure was necessary to pro- 
vide continuity over the interval between gage block 
lengths. There is too much uncertainty introduced 
by measuring an appreciable pressure difference 
directly by the diaphragm gage. However, the charac- 
teristics of the manometer cells are well known, so 
that an accurately known change of capacitance can 
be accurately related to a change of meniscus height. 
The bridge is balanced at "zero level" by adjustment 
of C Var , a 3-lead General Radio variable capacitor, Type 
1422-CE, 4 used in the range 0.005 to 0.11 pF, as shown 




4 Certain commercial materials and instruments are identified in this paper in order to 
adequately specify the apparatus and experimental procedure. Such identification does not 
imply recommendation or endorsement by the National Bureau of Standards. 



Figure 2. Interference pattern of gage block. 

One fringe equals about 0.25 /Am. 

in figure 3e. This capacitor has an imprecision of 0.0001 
pF and was calibrated at the values used with an un- 
certainty of the reference estimated at ±0.00003 pF. 
The uncertainties in the effective areas of the capac- 
itance plates ( ± 0.03% at 99% confidence) limit the 
attainable accuracy in calculating the change of menis- 
cus height from a change of C Va r- However, so long as 
the height of the mercury column differs from the 
gage block length by no more than 17 jLim, the contri- 
bution to the error in the measured vapor pressure 
from this cause will not exceed 1 ppm. 

2.2. Modified Triple Point Cell 

The vapor pressure apparatus, shown in figure 1, 
provided for circulation of the water. It was pumped 
from the pressure cell, PC, by a bubble pump, B into 
a reservoir, R. The reservoir was part of a still with a 
heater, H, to evaporate water into a condenser, C. The 
condenser water reentered PC at the top through a 
trap. The cell was 38 cm deep, with a volume of about 
500 cm 3 . A thermometer well, W, 32 cm deep, was a 
close fit on the 7 mm diameter case of a platinum 
resistance thermometer. Seven bulges, 20 mm across, 
supported the ice mantle when the water was pumped 
out of the cell. The bubble pump was a vertical tube of 
5 mm bore with a piece of perforated PTFE (polytetra- 
fluoroethylene) , PT, inserted in the lower half. A 50- 
watt strip heater was taped to the tube and thermally 
insulated from the bath. The reservoir was 50 mm in 
diameter and had a volume of about 700 cm 3 - It had 
a reentrant tube HP, which extended to the top of the 
boiler. The tube contained C2CI3F3 and functioned as 
a heat pipe, transferring heat from the heater, H, to the 
water surface. It was filled through the valve, V. There 
were three tubulations: one from the top of the cell to 
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Figure 3. Configurations of circuits used to measure the pressure. 



valve 4 was for the pressure measurements, and two 2.3. Pressure Transducer 

others, P, one from the side of the cell and one from 

the condenser, were for pumping from the gas phase. A special commercial pressure gage with a metal 

The tubulations were sealed by PTFE gaskets to stain- diaphragm, the position of which was sensed on one 

less steel high-vacuum valves at the top of the assembly, side by the capacitance between it and an electrode, 
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was used as a null device to detect equality of pressure 
between the water vapor on the one side and the 
helium gas transmitting the measured pressure on the 
other side. The null reading could be determined before 
and after measurements by opening the bypass, valve 2 
in figure 1, and adequately evacuating both sides. 
Great care had to be taken in the use of the instrument 
to achieve the needed stability, both electronic and 
physical. 

3. Filling of Vapor Pressure Apparatus 

The vapor pressure apparatus was cleaned and filled 
by the same techniques used in the preparation of 
conventional triple-point cells. The glass was cleaned 
with chromic acid, rinsed, lightly etched by hydro- 
fluoric acid, and leached by steaming. Two vapor 
pressure cells were filled in the course of a production 
run of triple point of water temperature standards, with 
water that was purified first by treatment with an ion 
exchanger and then quadruply distilled, with chemical 
treatment for removal of organic material between 
each distillation. 

Samples of the purified water were taken for isotopic 
analysis at the time of filling and also later during 
operation. Up to the present, there is no measurement 
of the vapor pressure of the triple point of naturally 
occurring water for which the maximum reported 
variation of isotopic composition would make a sig- 
nificant difference. However, the high accuracy of the 
present measurements being reported in this paper 
require such a specification to avoid added uncertainty. 

4. Procedure 

The operation of the equipment consisted of prepar- 
ing the manometer and the vapor pressure apparatus 
for measurements, determining the null of the evacu- 
ated diaphragm gage with the bypass valve open, and 
then, after closing the bypass valve, back filling the 
diaphragm with helium on the side of the capacitor 
electrode and with water vapor on the other. The vapor 
pressure apparatus was then vacuum pumped for 30 s 
in each of three successive 2 min periods before any 
actual measurements were made. We shall expand on 
these preliminaries and then describe the measure- 
ments themselves. 

Th » manometer is "zeroed" with a fixed mercury- 
to-capacitance-plate separation in the upper cell and 
for the corresponding separations in the lower cells 
when all three mercury menisci are on the same level. 
Under these conditions, the "zero level" capacitance 
necessary to balance the capacitance bridge (see 
fig. 3e) is established. The measurement of a vapor 
pressure consists of generating and measuring in the 
manometer the equivalent pressure of helium, as 
evidenced by a zero pressure difference at the dia- 
phragm pressure transducer. To do this, a gage block 
of appropriate height is inserted between the upper 
cell and its pedestal, and the pressure of helium in the 
lower cells and the volume of mercury in the manom- 



eter are adjusted until all mercury menisci are at a 
level relative to the capacitance plates to produce 
the same capacitances as for zero level. In these 
measurements, for the first time, we have modified 
this procedure in order to interpolate between the 
gage block values. We adjusted the manometer pres- 
sure until the diaphragm was balanced, and then we 
balanced the manometer capacitance bridge by 
adjusting C Var - 

The vapor pressure apparatus when not in use is 
operated in a "Standby" mode, in which the water is 
continuously circulated by pumping it out of the triple 
point cell into the still and distilling it back into the 
triple point cell. Occasionally, the vapor is vacuum 
pumped, the purpose being to remove gases released 
from solution, The operating parameters, and even 
some details of design, of this apparatus were deter- 
mined by long experimentation. The heat pipe in the 
condenser reservoir was necessary to avoid over- 
heating the water near the bottom of the reservoir. 
With the apparatus in a 0°C thermostat and a heater 
power of 25 watts, the water temperature in the reser- 
voir was about 25 °C and nearly uniform up to 1 cm 
below the surface. Distillation occurred at the rate of 
about 0.01 g/s. The bubble pump, operated with a 
heater input of 15 W, initially transferred a mixture 
of water and vapor from the full triple point cell to 
the reservoir at a rate of about 10 g/min, but the rate 
dropped as the water level dropped. When the cell 
was nearly empty, the pump acted as a still with a 
transfer rate of a few g/h. 

To prepare for measurements, the apparatus was 
placed in a bath, thermostated by circulating liquid 
close to the triple point temperature and constant 
within 0.2 mK. The vapor pressure cell was filled by 
distillation from the condenser (operated at 20 W this 
required about 24 h). Bath fluid was circulated through 
the condenser, and was pumped over the top of the 
cell and also onto the stem of a platinum resistance 
thermometer inserted into the well of the vapor pres- 
sure cell. (When water was being pumped out of the 
cell, a zone at the top of the reservoir was cooled by 
pumping bath fluid over it also). 

The apparatus was removed from the bath and a 
mantle of ice was frozen around the thermometer 
well using powdered solid C0 2 as the coolant. The 
completed mantle was from 5 mm to 10 mm thick. 

The apparatus was then returned to the bath, which 
was regulated at a temperature not more than 1 mK 
below the triple point temperature and the bubble 
pump was started. By the next morning, the cell was 
nearly empty whereupon the reservoir heater was 
turned on. It distilled water into the pressure cell at 
about 0.01 g/s (with the bubble pump left on, measure- 
ments could be conducted for several hours before 
enough water accumulated in the cell to touch the ice 
mantle). 

In the meantime, both sides of the diaphragm gage 
had been evacuated overnight with an ion pump. A 
"zero" reading was obtained with the gage isolated and 
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the bypass valve open. Then the bypass valve was 
closed, the electrode side of the diaphragm cell was 
filled with helium, and the other with water vapor 
(with the differential pressure across the gage kept 
within ±1 mm Hg). Finally, adjustments of the mano- 
meter variable capacitor were made in order to balance 
its bridge when the manometer pressure was such that 
the reading on the diaphragm gage was very near the 
"zero" value. 

Next, to sweep out gaseous impurities, the water 
vapor was pumped from the pressure cell, at a rate 
of approximately 100 cm 3 /s (0.6 cm 3 S.T.P./s). When 
the apparatus was in regular use, an initial set of three 
pumpings, each pumping lasting 30 s over a period of 
6 min was sufficient to insure no effect after recovery 
to a steady state would be observed with further 
pumping. 

Chart records were made of the diaphragm gage out- 
put for its zero, and immediately following the last 
pumpout, readings of the output of the diaphragm gage 
were recorded. Because the stirring of the thermostat 
had to be stopped to avoid shaking the diaphragm, 
the recordings were made for 20 s periods every two 
minutes. After each set of 3 recordings, the water 
vapor was again pumped for 30 s. After 4 such sets, the 
apparatus was closed off, the diaphragm bypass valve 
opened, and the diaphragm pumped out for 20 min 
when a new "zero" was recorded. During each mea- 
surement, the temperature of the diaphragm was 
sensed by reading the e.m.f. of a copper-constantan 
thermocouple referred to ice. 

5. Measurements 

Besides general experimentation to establish effec- 
tive operating techniques, there were three periods 
of several weeks during which actual pressure determi- 
nations were made. The first period was purely instruc- 
tive and resulted in vigorous efforts to improve various 
aspects of the equipment, in particular the diaphragm 
gage. The results of the second period reflected sub- 
stantial improvement in the accuracy of the measure- 
ments, and were reported at the 1974 meeting of the 
International Association for the Properties of Steam 
in Giens, France [8]. The results of the third period 
are decidedly more accurate than the second, more 
because of further improvement in the diaphragm gage 
than anything else. However, numerous possible effects 
were investigated in addition, and these results will 
be described here before presenting the final sets of 
measurements. 

Three standard calibrated platinum resistance 
thermometers were used to measure temperature, 
one in the well of the vapor pressure cell, the second 
to measure the temperature of the thermostat bath 
close to the vapor pressure cell, and the third as a 
regulator sensor. They were calibrated in a triple point 
of water cell each day; it is believed the values of the 
measured temberatures were uncertain by no more 
than 30 /*K. The temperature registered by the 



thermometer in the well of the vapor pressure cell 
could be used as an indication of the state of the 
mantle. When the ice mantle was dry or nearly so, 
pumping resulted in a rapid reduction in the pressure, 
by about 100 Pa, and in the temperature, by about 
0.1 K. When an adequate flow of water was maintained 
over the ice during pumpout, the drop in pressure was 
no more than 0.5 Pa and in temperature no more than 
0.1 mK. 

It has also been observed during a run earlier than 
any being reported that when air in an amount > 100 
mm Hg in pressure had entered the cell, the tempera- 
ture at the mantle was depressed by several tenths of 
a mK and did not recover fully before a new mantle 
was frozen. 

The effects on the observed vapor pressure resulting 
from variation of the bath temperature from the 
temperature of the triple point were studied. For 
temperatures in the range from 1 mK below the 
temperature of the triple point to 200 mK above it, 
scarcely any variation of the vapor pressure was 
observed. This indicated that the temperature sur- 
rounding the apparatus was not very critical, and fur- 
thermore, there was probably no undesirable effect 
from the warmer parts of the tubulations which ex- 
tended out of the bath. In the final measurements, the 
bath was thermostated at, or not more than 1 mK 
above, the temperature of the triple point of water. 

The diaphragm unit itself was not thermostated, 
but was thermally lagged with a thick layer of insula- 
tion. Its temperature was measured at frequent inter- 
vals, and in its final physical and electronic configura- 
tion, a consistent relationship between the copper- 
constantan e.m.f. 's and the zero readings existed. A 
straight line was fitted to the data by the method of 
least squares with the result that the zero can be 
calculated from the e.m.f. by the following equation: 
7=893.627-0.748503 V where y is the chart reading 
for the diaphragm zero and V is the thermocouple 
e.m.f. in \jS \ The residual standard deviation was 
1 chart division, and the standard deviation of a pre- 
dicted point was about 0.5 chart division. This equation 
was used to calculate the chart zero pertinent to any 
given chart pressure reading, from thermocouple 
e.m.f.'s observed at the same time. 

No effect dependent on refreezing the mantle was 
observed, and in fact, so long as crystals of ice re- 
mained, it appeared that the pressure was established, 
except for the change in the "pressure head", i.e., 
the pressure due to the weight of the column of vapor 
above the effective location of the line at which triple- 
point conditions exist. 5 

The values given in this paper were measured on 
four successive days and consist of the following five 

a . 11 

sets : 

1. Set 1: Ten values of diaphragm gage readings in 
3 groups, with zeroes before and after the 
readings. 



5 The calculated variation of the temperature of the ice mantle over its length because of 
the fluid head was less than 1 /iK. 
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2. Set 2: Nine values of diaphragm gage readings in 
3 groups, with zeroes before and after the 
readings. 

3. Set 3: Eight values of diaphragm gage readings in 
3 groups, with zeroes before and after the 
readings. 

4. Set 4: Twelve values of diaphragm gage readings 
in 4 groups, with zeroes before and after the 
readings. 

5. Set 5: Twenty-five values of diaphragm gage read- 
ings in 8 groups, with zeroes before and after 
each 4 groups. 

A group comprises the values measured between 
pumpouts. After the initial three 30 s pumpouts, the 
apparatus was pumped about every 10 min. At the con- 
clusion of all the measurements, the value of the 
capacitance to balance the capacitance bridge for the 
manometer at zero level was redetermined. 

6. Equations and Calculations 

Because of the effect of the pressure head, there is 
only a surface, a horizontal plane except for temper- 
ature perturbations, at which the triple point pressure 
can exist. The pressure of the water vapor at the line 
formed by the intersection of the surface with the 
mantle is given by the following: 



*tv — * \ 



8P V 



;+8P< 



+ 8P. 



Thermp 



+ oP Dia + 8P v 



ap Hd- 



^Man is given by eq (27) of ref. [7], with tempera- 
tures for zero level equal to 20 °C and for helium 
filling gas. The vapor pressure of mercury and its 
temperature derivative have been reevaluated at 
20 °C as 0.171 Pa and 0.0147 Pa K 1 , respectively, 
and the density at zero pressure is evaluated for 
* 6 8 = 20°C (rather than J 48 ), as 13,545.82 Kg/m 3 
[9]. Symbolically, 

^Man = - otpg (^uc & uc - h LC Si LC ) 



+ I iji JHgVap (o^uc 

+ pgh 2 (1 - a [8t Hli ]) + P 20 (Hg) 



20) 



+ pgh 2 (2.27 X 10" 6 h 2 



2.5 X lO" 7 

- 4.07 X lO" 3 Mho) , 



where 



• l ( dV \ 



for mercury, p is the density of 



mercury under zero pressure at 20 °C, g is the ac- 
celeration due to gravity, the dt's are differences of 
temperature from 20 °C of the upper cell (UC), the 
average of the lower cells (LC), and the mercury 
arm (Hg). The height of mercury in all the cells is 
substantiallv the same_ and is symbolized by A uc 
for the upper cell or by h LC as the average for the lower 
cells. The height of the mercury column is h 2 , and the 
pressure head is produced by the height of the gas 
column between the menisci of the lower cells and the 
center of the diaphragm (which was mounted ver- 
tically), designated ho. The pressure head depends 
upon the molecular weight M. The numbers are con- 
sistent for values of the symbolic quantities in S.I. 
"Manometer readings," comprising the observa- 
tions of the resistance of a capsule platinum resistance, 
thermometer located in a thermocouple reference 
block and the e.m.f.'s of all 4 sets of differential 
thermocouples for the upper cell, the mercury arm, 
the left lower cell and the right lower cell, were made 
once or twice for each set of vapor pressure measure- 
ments. The temperatures of the vertical components 
of the mercury lines were calculated from the thermo- 
couple e.m.f.'s as a difference from the temperature 
of the reference block. There were twelve copper- 
constantan thermocouples in series fastened to the 
upper cell, twenty-five to the mercury arm, ten to the 
left lower cell and twelve to the right lower cell. The 
observed e.m.f.'s were converted to temperature 
differences by a "handbook" sensitivity of 40.5 /iV 
°C~\ per junction at 20 °C, viz., 486, 1013, 405 and 
486 /xV °C _1 , respectively. The platinum resistance 
thermometer qualified as a standard instrument for 
the IPTS, with a reading in the triple point cell of 
/? PRT //? std = 0.25557275, and calibrated values of 



and 



A = 3.9855387 X 10" 3 °C" 



£ = - 5.8755669 X 10 ~ 7 °C- 2 



The values of temperatures were calculated by the 
program BRIDGE, given in appendix I. The following 
results (table 1) were obtained when the experimental 
data were calculated by the program PRSURE, 
given in appendix II. 

The manometer zero was established from 10 de- 
terminations made over a period of 10 days. The values 



Table 1. Manometer pressure at the diaphragm on the manometer side 
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were adjusted to that which would have been ob- 
served if all parts of the manometer were at a uniform 
temperature of 20 °C. The equation is derived from 
eq (27) of ref. 7, and as programmed in Basic is given 
in appendix III. The results in Pa were converted to 



nm of Hg and in turn the value of the variable ca- 
pacitor was adjusted on the basis that a change of 
1 nm is equivalent to 1 X 10" 5 pF. The values meas- 
ured and the calculated results were as follows 
(see table 2): 









Table 2. 


Observations of manometer zero 








PRT 














Corr 


Date 


Capsule 
(ft) 


0*V) 


v LCL 

(/*V) 


V LC K 

(/xV) 


£- Var 

(pF) 


CalcP 
(Pa) 


Corr 
(nm) 


Cyar 

(pF) 


11/27 
11:21 


27.583362 


-16 


-5 


-21 


0.02323 


4.04 X lO" 3 


+ 30 


0.02353 


1:45 


27.581889 


-17 


-6 


-18 


.02362 


7.15 X lO" 3 


+ 54 


.02416 


2:55 


27.581750 


-3 


3 


-8 


.02428 


1.29 X 10- 3 


+ 10 


.02438 


4:15 


27.581966 


-9 





-12 


.02409 


4.88 x lO" 3 


+ 37 


.02446 


5:17 


27.581802 


-14 


-3 


-16 


.02391 


6.83 x 10- 3 


+ 51 


.02442 


12/3 
10:10 


27.593201 


-17 


-6 


-24 


.02356 


2.24 X lO- 3 


+ 17 


.02373 


2:03 


27.591870 


-18 


-6 


-21 


.02357 


6.33 x 10- 3 


+ 48 


.02405 


12/4 
9:25 


27.587772 


-17 


-5 


-20 


.02343 


6.49 X lO 3 


+ 49 


.02392 


12/5 
9:50 


27.584449 


-17 


-5 


-19 


.02346 


7.31 x 10- 3 


+ 55 


.02401 


12/6 
11:41 


27.582260 


-15 


-5 


-18 


.02333 


4.87 X lO- 3 


+ 37 


.02370 




Ave 0.02404 
















Calib. Cc 


>rr+ 0.00031 




Net 0.02435 



The estimate of the standard deviation of the mean is S = 0.00012 pF. 



S^wring is the net correction for imperfectly joined 
gaging surfaces between the wringing boss on the 
bottom of each cell and the pedestal or the gage block. 
The leak rate through the gap formed when two 
gaging surfaces meet on one edge with a dihedral 
angle </> (in /x rad) was studied [10] and is given for the 
particular blocks and volume of vacuum system as 

dp/dt = 0.42 (f) 2 + 0.022 3 /xm Hg/min. 



arately. It was calculated that the volume of each cell 
plus thermocouple gage is related to the total volume 
from which the original equation was derived as 

VjjcJV Tot = 0.579 

V L cnlV Tot = 0.620 



^uc/^t 



0.523. 



The vacuum leak rates were measured, both The rates were remeasured after every adjustment of 

during the determination of zero level and several any cell position. The net pressure effect is given in 

times during the course of the vapor pressure measure- the following: 
ments. The gaging surfaces are square, of a width 



W = 2.41 cm on a side. The change in height between 
the measurement of zero level and the pressure meas- 
urement was then 

8h = [</>(UC) - </>o(UC) - {<MLCL) + 0(LCR) 

- </>o(LCL) - <MLCR)}/2] X W\l 

and, with dh in nm, the pressure effect in pascals is 



Table 3. Pressure effect from imperfect gaging joints 



8P, 



Wring ' 



dh x 1.32763 x 10- 4 Pa. 



The leak rates were measured for each cell sep- 



Set 


dp/dt (/xm Hg/min) 


$P Wring 


LCL 


LCR 


UC 


(Pa) 


Zero Level 
1,2 and 3 
4 and 5 


15 
16 
16 


39 
43 
40 


5 
23 

22 


0.00373 
0.00357 
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S^cvar is derived from the change of mercury level 
requiring a change of the variable capacitor to main- 
tain the manometer bridge balance. The change of 
level can be accurately determined because of the 
coaxial switching arrangement that is part of the 
manometer instrumentation. The configurations used 
are shown in figure 3. The manometer transformer has 
two taps on the right-hand side, and one on the left. 
The tap /?i has 3/4 the voltage of L, and R> has the 
same voltage as L. The variable capacitor may be 
switched to any tap. 



There is also a switch which changes ground from 
the center tap to the other side of the detector. There 
are 3 standard capacitors, labeled C u C Ut and C m . 
The values of C, and C„ are 8.2046 ±0.0008 pF, and 
are equal within 1 ppm. The value of C m is 3.0503 
pF, so that it can be used to measure the lower cell 
capacitances. The capacitances of the various parts 
of the detection system and the manometer were 
found from the following sets of measurements at 
zero level: 



Table 4. Measurement scheme to determine manometer capacitance characteristics 





Switch 




Center 






Figure 


Pos. 


Tap 


Tap 


^ Var 


Capacitance Equation 


3a 


I 


Ri 


Gnd 


Out 


C j = 6n 


3b 


II 


Ri 


Gnd 


Out 


C a — C i 


3c 


I 


R, 


UnGnd 


0.04313 pF 
Corr. 0.00055 pF 


C f "T C i — C jj i U Var 

C T = 0.04368 pF 


3d 


III 


Ri 


UnGnd 


Out 


C x -+• C uc = c. ii 
Cue = 8.16092 pF 


3e 


IV 


Rx 


UnGnd 


0.02373 pF 
Corr. + 0.00006 pF 


C T + SGlc = 3/4 (C vc + C Var ) 
ICuc = 6.13853 pF - 0.04368 pF 




.02379 pF 


2C IjC = 6.09485 pF 


3f 


VI 


Ri 


UnGnd 


0.01215 pF 
Corr. + .00012 pF 

0.01227 pF 


C T + C ixr + C var = C i,cl 

C lC R=C in -0.05595 pF 
C LCI = 3.07540 pF 
C i CR = 3.01945 pF 



The mercury cell capacitors comprised guarded 
circular plates facing a grounded surface that was 
essentially infinite. The capacitance can then be ex- 
pressed as C= 1.11267 A/(4 tt S) pF, where A is 
the effective area of the plate (in cm 2 ), and S is the 
separation between the plate and the mercury (in 
cm). The dimensions of the cells are accurately 
known; the following separations at zero level can 
then be calculated (table 5): 

Table 5. Characteristics of manometer cells 



was measured with the variable capacitor in the L arm, 
so that C T + C LC + C Var — 3/4 C uc and the change of 
capacitance of the lower cells from the balance at 
zero level is 0.75 X 0.02435 pF plus the reading of 
C V ar f° r tne Set. This change of capacitance corre- 
sponds to a displacement, d 9 in the separation of the 
mercury surfaces from the plates in the lower cells, 
which is 



AXC LR = 



1.11267 





Cap. a 


Eff. Diam. 


So 


Cell 


(pF) 


(cm) 


(/im) 


UC 


8.16092 


3.07838 


807.517 


LCL 


3.07540 


1.88488 


803.364 


LCR 


3.01945 


1.88509 


818.432 



A ( I I ^ 

LCL \S (LCL)+d So(LCL)/ 
+ Alcr \ So(LCR) + d ~ So(LCR) j 



a The number of digits is consistent with the imprecision. The total 
uncertainty is nearly 1 part in 10 4 , but the net changes are known 
within 3 ppm of the total figure. 

In operation, the mercury level is maintained so that 
Cue — 8.16092 pF. The value of the variable capacitor 
to balance the bridge at zero level, the "manometer 
zero," corrected for differences of mercury tempera- 
ture from 20 °C, was 0.02435 pF. The vapor pressure 



The pressure in pascals is related to the displacement 
in jLtm as 

Pc V ar= 0.132763 X d Pa. 

The values obtained are given in table 6. 

Table 6. Pressure increment determined by capacitor setting 



Set 


Cvar 

( P F) 


AC 

(pF) 


d 
, (/Am) 


oP CVar 

(Pa) 


1 
2 and 3 
4 and 5 


0.02644 
.02531 
.025306 


0.044703 
.043573 
.043569 


5.9911 
5.8380 
5.8376 


0.79532 
.77507 
.77502 
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8P Thermp , the effect of thermomolecular pressure, 
is expected to be significant only in the pressure trans- 
mitting tube of the manometer. This tube, which is 
0.1512 cm in radius, ran from the manometer vault, 
at a temperature very close to 20 °C, to the valves that 
were physically adjacent to, and assumed to be at 
the temperature of, the diaphragm unit. 

The value of the effect was derived from the equa- 
tion of Weber, Keesom and Schmidt [11] (the program 
in Basic for their equation is given in appendix IV). 
We measured the thermomolecular pressure dif- 
ference of helium in a tube of 0.0412 cm radius and 
found the values to be less than the values from 
Weber's equation by a factor of 0.76. The functional 
relationship in the Weber equation for the radius of 
the tube, i.e., that pAp oc 1/r 2 , is assumed to be correct. 
The difference between theory and experiment is 
thought to occur because of the inadequacy of other 
assumptions. Therefore, the same factor of 0.76 was 
used in the calculation because it should still be 
applicable for the larger tube. (See table 7.) 

Table 7. Thermomolecular pressure correction 







Ap (Weber) 


0.76 Ap 


Set 


(°C) 


(Pa) 


(Pa) 


1 


25.05 


0.02122 


0.01613 


2 


25.13 


.02155 


.01638 


3 


25.13 


.02155 


.01638 


4 


25.03 


.02105 


.01600 


5 


24.82 


.02021 


.01536 



SPdia> the pressure difference at the diaphragm, 
was derived from sets of charts readings. Following 
the stated operating procedures, we recorded the dia- 
phragm zero and the copper-constantan thermocouple 
reading, then pumped three times, after which three 
sets of measurements, each set consisting of a thermo- 
couple reading and a diaphragm reading for the pres- 
sure difference, were recorded. The diaphragm 
output was recorded for 20 s every 2 min. The thermo- 
couple was read again while the water vapor was 
repumped from the vapor pressure cell. Three more 
thermocouple and pressure readings were taken for 
20 s at 2 min intervals, and then the diaphragm was 
evacuated for 20 min and a diaphragm zero recorded. 
It was experimentally demonstrated that this period 
of evacuation of the diaphragm was long enough to 
assure reliable zero values. 

Within the imprecision of the instrument, the 
diaphragm readings were correct on the 10 /mm Hg 
scale, so that the chart measurements could be cali- 
brated by the diaphragm readings. In order that 
the uncertainties would be small, the manometer 
was adjusted so that the difference of the diaphragm 
from zero was small. The net deflections in chart 
divisions in Set 1 were — 2, — 2.2, — 1.7, —6.6, — 4.7, 

- 1.7, — 5.5, — 3.7, —4.7, and — 5.4, for an average of 

- 3.82 (S = 1.65). For set 2, they were - 0.2, 0.8, 

- 0.8, 0.8, 2.8, 2.5, 2.5, 2.5 and 1.5, for an average of 



1.38 (S = 1.31). For set 3, they were 1.1, 0.1, 0.9, 1.5, 
0.5, - 1.5, - 1.2 and - 0.8, for an average of- 0.12 
(S = 1.04). For set 4 they were 1.1,2.2,-0.1,-2.1, 
1.2, - 1.7, 2.8, 4.3, 0.3, 1.3, - 0.1, 0.6, for an average 
of 0.82 (S = 1.79). For set 5 they were 2.5, 1, 3.4, 0, 
2.5, 1.5, - 0.5, 1.5, 1.5, - 0.5, 2.2, 1.2, 3.1, 2.1, 
0.2, - 1.8, - 1.3, 0.2, 1.1, 1.1, 0.6, 2.6, - 1.9, 0.2 
and 1.6, for an average of 0.89 (S= 1.39). 
The chart calibration gave a sensitivity of 0.00536 
Pa/div so that we have for the five sets (with S denoting 
the estimate of the standard deviation for a single 
measurement and S m the estimate of the standard 
deviation of the mean) (see table 8): 

Table 8. Measured pressure differences at the diaphragm 





Ave Defl 


8Pdi\ 


5 


S 


Sm 


Om 


Set 


(div) 


(Pa) 


(div) 


(Pa) 


(div) 


(Pa) 


1 


-3.82 


-0.0205 


1.65 


0.0089 


0.52 


0.0028 


2 


1.38 


.0074 


1.31 


.0070 


.43 


.0023 


3 


-0.12 


-.0006 


1.04 


.0056 


.35 


.0019 


4 


.82 


.0044 


1.79 


.0091 


.52 


.0028 


5 


.89 


.0048 


1.39 


.0075 


.28 


.0015 



8 ^VapHd? tne pressure head produced by the column 
of water vapor, is calculated from the estimated verti- 
cal distance, h, from the midpoint of the diaphragm to 
the line at the mantle where the triple point is realized. 
We believe this line to be near the top of the mantle, 
but its position is not entirely stable, both because the 
flow of water may vary and because the mantle melts 
and recedes from the top. The pressure is calculated as 



8P, 



VapHd 



^Man g M 

RK 



Ti 7W' 



where g is the acceleration due to gravity, M is the 
molecular weight of the vapor, R is the molar gas con- 
stant, K is the conversion factor between units of pres- 
sure, hi is 0.244 m running between the diaphragm at 
25 °C and the bath at 0°C, and h 2 is the distance, 0.07 
m, in the bath at °C. Then 



°AapHd = 23.31 X 10- 6 P* 



0.01424 Pa. 



The average temperature, 7\, was weighted on the high 
side. The value is bounded by the extremes, of course, 
and within those limits any reasonable choice would 
hardly make any significant difference in the final 
pressure. 



7. Results 

The results are summarized by totalling the elements 
of pressure for each set (all values in pascals): 



514 



Table 9. Total values of the vapor pressure of water at its triple point 



Set 


'Man 


S^Wrinn 


8'Cvar 


°*Thermp 


S^DIA 


0* VapHd 


Ptv 




(Pa) 


(Pa) 


(Pa) 


(Pa) 


(Pa) 


(Pa) 


(Pa) 


1 


610.84306 


0.00373 


0.79532 


0.01613 


- 0.0205 


0.01424 


611.6520 


2 


610.84504 


.00373 


.77507 


.01638 


.0074 


.01424 


611.6619 


3 


610.84572 


.00373 


.77507 


.01638 


-.0006 


.01424 


611.6545 


4 


610.84403 


.00357 


.77502 


.01600 


.0044 


.01424 


611.6573 


5 


610.84662 


.00357 


.77502 


.01536 


.0048 


.01424 


611.6596 



The average of the 5 sets is 611.6571 Pa (S = 3.94 mPa). The estimate of the standard deviation of the mean isS m = 1.76 mPa. 



8. Isotopic Composition 



The determination of the isotopic content of samples 
from the vapor pressure cell relative to SMOW 
(standard mean ocean water) depended upon two steps. 
First, the isotopic content of cell water was deter- 
mined relative to NBS-1 standard water [12], the 
composition of which had, as a second step, been 
related to SMOW by the work of Craig [13]. The 



DjH absolute abundance ratio of SMOW has been 
determined very precisely by Hagemann, Nief and 
Roth as 155.5 ppm [15]. The value for the absolute 
abundance of 18 in SMOW is based on Craig's 
measurements [13, 14] and as variously interpreted 
may be taken to be 18 /0 16 = 1993.4 X 10~ 6 to 1995 X 
10~ 6 , of which the latter is used here. The composi- 
tions of seven samples of the water used in the vapor 
pressure cell, both relative to SMOW and absolute, 
are as follows: 



Table 10. Isotopic composition of water in the vapor pressure cell 





Sample 

No. 


Description 


D/H% 


D(ppm) 


018/016 % 


,K (ppm) 


1st cell 
2nd cell 


1 
2 
3 

4 
5 
6 

7 


raw (tap) water 
Beginning of cell filling 
End of cell filling 

Beginning of cell filling 
raw (tap) water 
End of cell filling 

Distilled from apparatus 2 weeks after 
experiment 


-4.0 
-4.9 
-5.6 

-4.6 
-5.5 
-5.1 
-6.9 


149 
148 

147 

148 
147 
148 
145 


-0.58 
-0.71 
-0.88 

-0.60 
-0.77 
-0.64 


1983 
1980 
1977 

1983 
1980 
1982 


Average 


147 


1982 















All measurements were made with the second cell. 
The difference between the vapor pressure of the cell 
water and SMOW because of the reduced compositions 
of D and 18 can be calculated from measurements of 
Majoub [16], who reported that the ratio of / > hdo/^h 2 o 
at 0°C is 0.8994, and of various authors [16, 17, 18], 
who reported ratios of Ph 2 o™IPh 2 o™ at °C that 
average about 0.9885. The difference of pressure may 
then be derived from Raoult's Law, where the pressure 
of the liquid of the sample at 0.01 °C is 1 ppm more 
than the corresponding metastable SMOW. The 
sample triple point temperature is about 40 fxK lower 
than the triple point of SMOW [19]. Given that dP/dt 
is about 44.4 Pa/K for water at 0.01 °C, the effect of 
the difference in temperature of the triple points is 
about 3 ppm. Thus the vapor pressure at the triple 
point of the sample, which corresponds closely to the 
analyses of typical "continental water," can be cal- 
culated to be about 2 ppm, or 1.2 m Pa, less than the 
vapor pressure at the triple point of SMOW. 



9. Discussion 

It is, of course, crucial to the reliability of any data 
that procedures be established to eliminate or satis- 
factorily minimize possible errors before any final 
measurements are made. We believe that the probable 
difficulties requiring special study were (1) the pos- 
sible variation of the vapor pressure because of a 
slow return to equilibrium following a pumpout, or 
with variation of ambient temperature, (2) the presence 
of contaminating gas in the system, and (3) the pos- 
sible errors in determining the diaphragm zero. 
Each of these problems received special study. 

As discussed in section 5, the pressure and tem- 
perature were only slightly perturbed by pumpout 
so long as all three phases were maintained about the 
mantle. Similarly, in the presence of the mantle, 
there was at most only a slight change in vapor pres- 
sure when the ambient temperature varied from 1 
mK below to 0.2 K above the triple point. Thus the 
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actual realization of the triple point must have con- 
tributed great stability to the vapor pressure; the 
broad expanses of exposed liquid and solid phases 
could be expected to fulfill the physical necessities 
for facilitating such an equilibrium. 

The procedure of pumping for three 30 s intervals 
prior to measurements, and subsequently for 30 s 
every 10 min was adopted because the measured 
pressure, with lesser amount or frequency of pumping, 
varied in a way to indicate the accumulation of meas- 
urable amounts of contaminating gases. The three 
recorded diaphragm pressure readings made in each 
interval between pumping showed no evidence of 
preferential drift that might indicate the rate of 
pressure buildup to be significant in that span of 
time. 

An accurate measurement of the diaphragm zero 
depended upon the evacuation of enough of the sorbed 
gases (particularly water vapor) that the remainder 
would not exert a significant pressure difference 
between the two sides of the diaphragm, with the 
by-pass valve open. It was observed experimentally 
that a 20 min evacuation following pressure measure- 
ment appeared to be sufficient. These measurements 
together with the thermocouple e.m.f.'s observed at 
the same time correlated with those values of zeroes 
after long pumping and their concomitant thermo- 
couple e.m.f.'s. This fact offered substantial confir- 
mation both that 20 min was an adequate pumpout 
and that the diaphragm zero could be accurately 
associated with the thermocouple readings. However, 
measurements made when the room temperature 
varied rapidly were not in agreement with the rest of 
the results, and were not used. The assumption was 
that the thermocouple, being on the outside of the 
diaphragm, could not have represented the appropri- 



ate temperature when there were substantial gradients. 
The reader may note that the corrections for 
C V ar appear to be inconsistent. This is because a 
different reference line for reading was used at dif- 
ferent times. 

10. Estimation of the Total Uncertainty 

The stated total uncertainty consists of the limits, at 
99 percent confidence level for the random errors; 
and the systematic errors, estimated conservatively 
enough to warrant about the same confidence level. 

The estimates of the values of the errors are ex- 
pressed initially as one standard deviation. It is useful 
to treat them in two groups: (A) Those sources of 
error which contribute to the imprecision calculated 
from table 9, and (B) those additional sources of error 
which have a fixed effect on the directly determined 
value. 

The term "systematic error" can be understood, and 
is used by some authors, to denote any error that in- 
troduces a constant bias into the observed results. 
For example, this is true of each of the items of Part 
B. However, in all but the last entry the errors were 
evaluated from well-defined imprecisions, and do not 
differ in their nature from the imprecisions in Part A. 
These imprecisions are combined in a sum of the 
squares addition with those of Part A. In a more re- 
stricted definition, we are referring to systematic 
errors as those which introduce a bias into the results 
but which are not well enough investigated to be evalu- 
ated by statistical techniques, or are not realized to be 
important. The limits of a systematic error are apt to 
be less well-defined than for a random error, although 
there must be some theoretical or experimental 
results upon which an estimate can be based. 



Table 11. Sources of error 

(A) Random errors Observable in the Final Values 
Estimate of the standard deviation in the pressure as a result of imprecision in thermocouple e.m.f.'s 
Estimate of the standard deviation in assignment of the diaphragm zero (1/2 div.) 
Estimate of the standard deviation of the pressure head because of uncertainty in the location of the 

triple phase line (within 2 cm at 99% confidence) 
Estimate of the standard deviation of the pressure determination (the weighted average of the standard 

deviations of the diaphragm readings for the 5 sets) 



(Pa) 



S = 0.O020 
S = 0.O027 

S = 0.0003 

S = 0.0019 



The square root of the sum of the squares of this group is 0.0039 Pa and agrees with the imprecision of the results in table 9 expressed as 
the estimate of a standard deviation of a single measurement. The estimate of the standard deviation of the mean for 5 sets is 0.0017 Pa. 



(B) Random Errors and Systematic Errors Not Observed as Imprecision in the Results 
Estimate of the standard deviation of the mean of the manometer zero (based on 10 determinations) 
Estimate of the standard deviation of the mean of the gage block calibrations (based on 20 determina- 
tions) 
Estimate of the standard deviation of the mean of calculating the effect of imperfect gage block joints 

(based on the standard deviation of the calibration) 
Estimate of the error from the variable capacitor calibration (based on comparisons with a capacitor 

calibrated with 10 ppm total uncertainty) 
Estimate of the error in calculating the pressure from the variable capacitor readings (based on the 

measurement uncertainties of the apparatus dimensions) 
Estimate of the error in the calculation of the thermomolecular pressure effect (the random error based 
on the standard deviation of the calibration) 



S m = 0.0016 

S m = 0.0003 

S m = 0.0005 

5 = 0.0002 

S = 0.0005 

Random S m = 0.0003 
Systematic 0.0003 



The sum of the random errors in Part B, combined by the square root of the sum of the squares, is 0.0018 Pa. There is also a systematic error 
of 0.0003 Pa. 
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We attributed an element of systematic error to the 
last entry in Part B, inasmuch as neither theory nor 
experiment for thermomolecular pressure effects is 
complete. Another important systematic error might 
be the effect of contaminating gases on the measured 
pressure. We believe our equipment and procedures 
were devised in such a way that we were able to prove 
that this error could not have been significant in our 
measurements. 

The combination of the various random errors for 
99 percent confidence limits requires that they be 
weighted differently according to the number of 
degrees of freedom. The random errors in Part A 
were applicable to the 5 final pressures of table 9. 
With 4 degrees of freedom, the standard deviations 
were multiplied by 4.6 as given in Student's Table 
[20], to compute the limits for 99 percent confidence. 



In Part B, the number of degrees of freedom is nine or 
more, so that a multiplier of 3.2 was used. The separate 
weighted standard deviations of Part A and Part B 
were combined by the sum of the squares, to give an 
estimated uncertainty of ±10 mPa at the 99 percent 
confidence limits. The total additional systematic 
errors are estimated to be relatively insignificant. 



The authors express their gratitude to James L. 
Cross for his highly competent efforts in cleaning and 
filling the vapor pressure apparatus, to John S. Beers 
for his elegant calibration of the 4.6 mm gage block, 
and to Irving Friedman for accurately characterizing 
the isotopic compositions of the water samples. 



Appendix I 
Fortran IV program for calculating t,n from PRT measurements 



100 * FOPTPAN IV PPOGPAM NAMED EFI DGE FOP CALCULATION OF 

110 * IP??-6S 1'ALVES FPOM PPT PEFI FRANCE MEAFVPEMENTS . 

121? D0TELE PPECISION PTC 5) * PC 5) * W# AC 10) s EC 10) , ?C 5>* T1#Z 

130 DIMENSION N( 5) 

140 CALL DEFINE C 1, 7HTCAL74* ) 

150 CALL DEFINE C2,4HPPT, ) 

160 ^EAD (U110)(A(J)iE(J)iJ=U10) 

170 110 F0PMATCD15.8*2X*D15.8> 

130 120 PEAD C2j 130>CS1*S2*CNCJ)#PT<J>*PCJ>>J«1#4»>) 

19 130 FOPMATC 2X* F6 . 2, Xj F6. 3j ( 4C /2Y, I i\» D 1 5 . 8, D 1 5 .8 ) ) ) 

200 VPITE C9# 130>CSUS2*CNCJ)#PTCJ>*PCJ>*Jsl.>4>> 

2 10 VPITE C9j 150) Sl> S2 

220 150 F0PMATC/UX> , STATE , F6.2*2X/F6.3^/) 

230 DO 240 J= \j 4 

240 V=PT( J)/CPCJ)-0.996D-03) 

250 I=N(J)-700 

260 TCJ)*C-ACI >+D?GPT(ACI )**2-4.*C 1 •-V)*EC I > > ) /C 2 

27 ? = e.45D-0 3*' r <J)*(TC J)/ 100.-1.)*CTCJ)/419.58-1 

280 T <J)=T(J)+? 

290 VPITE C9*230)NCJ)*TCJ> 

300 230 FOPMATC 'THEPM. NO • ' I 3# 3X 

3 10 240 CONTINUE 

320 T l = CTC2) ♦ T (3) + TC4>>/3 
330 VPITE C9* 270)T1 

340 270 FOPMATC 'AVE. TCINT. 1968) 
350 GO TO 120 
360 END 



*ECI ) ) 
)*CT(J)/630. 



74- 1 . ) 



•TCINT. 1968) = ' D17.10) 



D17. 10,/S/) 



TCAL7I4 



0. 


39849575D- 


•02 


-0 


.5S762415D- 


■06 


0. 


39850801D- 


02 


-0 


,587721 19D- 


•06 


0. 


39848143D- 


02 


-0 


.58762135D- 


06 


0. 


398499 17D- 


■02 


-0 


► 58769654D- 


■06 


0. 


39849663D- 


02 


-0« 


.58756733D- 


•06 


P. 


39S47384D- 


02 


-0« 


►58757724D- 


■06 


0. 


39348898D- 


02 


-0. 


.58760947D- 


■06 


0. 


39850236D- 


2 


-0, 


58773636D- 


•06 


0. 


39838174D- 


02 


-0. 


.58757363D- 


•06 


0. 


39855337D- 


■02 


-0 


.58755669D- 


•06 
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Appendix II 

Basic Program for Calculating Pressure From Manometer Readings 

100 REM BASIC PROGRAM NAMED "PRSURE" 

110 REM TO CALCULATE EQ. 27 FOR MANOMETER. 

120 REM UNITS ARE SI* RESULTS IN PASCALS. 

130 REAO R0*G1*M*H0 

140 DATA 1.354582E 04*9 .801 022* 0. 40026E-2* 5. 1 76 

150 DATA 486*1013/405*486 

160 FOR I «1T04 

170 READ XU> 

180 NEXT I 

190 READ S1*?2*T0*H2 

200 FOR 1=1 TO 4 

210 READ E<I) 

220 NEXT I 

230 FOR I=1T04 

240 LET T(I)=E(I)/X<I) 

250 LET T<I)»T<I)+T0 

260 NEXT I . 

270 Tl*<T<3>+TC4>>/2 

280 P»0.8079*(T1-T( 1>> 

290 P=P+0.0147*(T< l>-20) 

300 P=P+R0*G1*H2*< 1-0. 17255E-3*<TC 2> -20) > 

310 LET P=P+R0*Gl*H2*<2.27E-6*H2-2.5E-7-4.07E-3*M*K0> 

315 P=V*» 171 

320 PRINT SUS2;P 

330 GO TO 190 

340 DATA 1.34* 1. 1* 19. 9389*4. 6000292E-3 

350 DATA -3* 2 1*4* -2 

360 DATA 1.34* 1.2* 19. 9395*4. 6000292E-3 

370 DATA -6* 23* 3* -4 

380 DATA 1.34* 1.3* 19 .9428* 4.6000292E-3 

390 DATA -5*25* 5* -3 

400 DATA 1.34* 1.4* 19. 9513*4. 6000292E-3 

410 DATA -6* 24* 4* -5 

420 DATA 1.34* 1.5* 19. 9530*4. 6000292E-3 

430 DATA -8* 2 1*4* -6 

540 END 



1*34 
1.34 
1.34 
1.34 
1.34 



1.1 

1.2 
1*3 
1.4 
1.5 



610.84306171 
610.84503579 
610.84572239 
610.84402799 
610.84661895 
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Appendix III 
Basic Program For Calculating Pressure Effects at Zero Level 

100 REM BASIC PROGRAM NAMED "ZPRES" 

110 REM TO CALCULATE EQ. 27 FOR MANOMETER ZERO* 

120 REM UNITS ARE SI* RESULTS IN PASCALS. 

130 READ R0/GUM/H0 

1/10 DATA 1.354584E 04*9 .80 1022, 0. 40026E-2* 5. 1 76 

150 DATA 486* 1013*405*485 

160 FOR I =1T04 

170 READ X(I) 

180 NEXT I 

190 U0 = 20 

200 READ S1*S2*T0 

210 FOR 1=1 TO 4 

220 READ EC I) 

230 NEXT I 

240 FOR I=1T04 

250 LET T(I)=ECI)/XCI) 

260 LET T(I>=TCI>+T0 

270 NEXT I 

280 Tl = (T(3>+T(4))/2 

290 P=0.8079*(T1-T< 1)) 

300 P=P+Q. 0147*CT< 1>-T1> 

310 PRINT si; S2;P 

320 GO TO 200 

330 DATA 1048*1.0*19.9567 

340 DATA -16/ 5* -5* -21 

350 DATA 1048*1.1*19.8094 

360 DATA -17* 7* -6*- 18 

370 DATA 1048*1.2*19.7855 

380 DATA -3* 19* 3* -8 

390 DATA 1048*1.3*19.8169 

400 DATA -9* 1 1*0*- 12 

410 DATA 1048*1.4*19.8 007 

420 DATA -14*8*-3*-16 

430 DATA 1048*1.5*20.0502 

440 DATA -17* 7* -6* -24 

450 DATA 1048*1.6*20.0369 

460 DATA -18* 6* -6* -21 

470 DATA 1048*1.7*19.9959 

480 DATA -17* 7* -5**20 

490 DATA 1048*1.8*19.9627 

500 DATA -17* 9* -5*- 19 

510 DATA 1048* 1.9* 19.9408 



520 DATA 


-15*7*-5*-18 






530 END 










1048 


1 


.404552E-02 




1048 


1.1 


•715303E-02 




1048 


1.2 


. 129202E-02 




1048 


1.3 


.487762E-02 




1048 


1.4 


•682625E-02 




1048 


1.5 


•224516E-02 




1048 


1.6 


•633001E-02 




1048 


1.7 


•649643E-02 




1048 


1.8 


•731340E-02 




1048 


1.9 


•48 6854E-02 
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Appendix IV 

Basic Program For Calculating Thermoinolecular Pressure 

100 PEM BASIC PPOGPAM "THEPMP"* DEPIVED FPOM VEEEP- SCHMI DT EQ. FOP 

110 PEM CALCULATING THE THEPMOMOLECULAP PPESSUPE DEVELOPED BY 

120 PEM HELIUM IN A TUEE OF PADIUS P PUNNING BETWEEN TEMPEPATUPES OF 

130 PEM Tl TO T2. WITS CGS, PPESSUPE IN DYNES/CMt2. 

140 PEAD L0iT0iKiiK2/N0/N 

150 PEAD P*M*P0 

160 DATA 17.85*273. 15* 1.36* 1.68* 1.86885E-04*0. 1866 

170 DATA . 1492* 1.25* 1.7846E- 10 

180 P=61 16.4 

190 PPINT "P =" P 

200 PEAD T2*T1 

210 Tl = Tl + 273. 15 

220 T2=T2+273. 15 

230 PPIN^ "Tl «" Tl*" T2 «" T2 

240 LET S=CT1/T2)**(M+ 1) 

250 LET S0=CT1/T0)**CN+ 1) 

260 LET L1=L0*S0 

270 LET E=4*K2*L1/P 

28 LET F=M- 1/C 2* <N+ 1 > > 

290 LET D=3*F*K 1/(2*3. 14159*K2) 

300 LET G=0.5*( (E/P) t2)*( 1+St 2)*( 1-D) 

310 LET H=2*E*C 1+St 2/ ( S+ 1) )/( 3*P> 

320 LET C=l-H+G 

330 LET U=( Tl**< 2+2*N)-T2**(2+2*N) ) / C T0** ( 2+ 2*N) ) 

340 LET V=6*Kl*(N0t 2)*U*C 

350 LET V=C l+N)*P0*Pt 2 

360 LET X =U/U 

370 PPINT "Plt2 -P2t2 = " X 

380 LET P3=X/C2*P> 

390 LET Y=P/(P+P3/2) 

400 LE^ P3=P3*Y 

410 PPINT "DELTA P = " P3 

420 PPINT 

430 GO TO 200 

440 DATA 19.94,25.05/19.94*25.13,19.96*25.03 

450 DATA 19.95/24.82 

460 STOP 

470 END 



7XBASIC THEPMP DE 

PUN 

P = 61 16.4 

Tl = 298.2 T2 = 29 3.09 

Plt2 -P2t 2 = 2595.4871 1 

DELTA P = .21217073 

Tl = 298.28 T 2 = 29 3.09 

Plr2 -P2t2 = 2636.4391679 

DELTA P = .2155 1834 

Tl = 298. 18 T2 = 293. 11 

Plt2 -P2t2 = 2575.2255458 

DELTA P = .2105 1446 

Tl = 297.97 ^2 ■ 29 3. 1 

Plt2 -P2t2 = 2472.7760644 

DELTA P = .202139 76 
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